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APPENDIX - RATE OF RISE OF GAS BUBBLES IN DISTILLED WATER 


NOTATION 


Projected area of body 

Drag coefficient (3 Fal 
3 y2 

Drag of body 

Coefficient of surface viscosity 


Acceleration due to gravity 


Length parameter of body 


gus 
Dimensionless parameter (- ot.) 

po 
Equivalent radius of bubble, i.e., radius of a sphere of equal volume 


zeae) 
pe 


Reynolds number (- 
2r.Up -) 


e 


Reynolds number inside fluid body (- 


Terminal velocity 


eee 


Weber number (- 

o 
Coefficient of dynamic viscosity of fluid medium 
Coefficient of dynamic viscosity of fluid inside fluid body 
Density of fluid medium 


Density of body 


Interfacial tension 


ABSTRACT 


In connection with other investigations at the David Taylor Model Basin, 
detailed information became necessary on the motion of air bubbles in variable 
pressure fields. Since no information on the subject was available, a fundamen- 
tal study of the motion of bubbles was undertaken. As an initial step, experi- 
ments were conducted to determine the drag and shape of single air bubbles 
rising freely in various liquids. 

The results of the experiments show that a complete description of the 
motion of air bubbles is not possible by use of dimensionless parameters con- 
taining the usual physical properties of the liquid (viscosity, surface tension, 
density). Three types of bubble shapes were observed in each liquid, namely 
spherical, ellipsoidal, and spherical cap. For a specific liquid, the shape of 
the bubble was a function of its volume. 

For tiny spherical bubbles, the drag coefficients coincide with those of 
corresponding rigid spheres. With increase in bubble size, a decrease in the 
drag as compared to that of rigid spheres occurs in some liquids. Thus, the 
drag curves of the spherical bubbles rising in various liquids fall between two 
limiting curves, namely the drag curve of rigid and fluid spheres, respectively. 
It was not possible to determine a criterion for the transition of the bubbles 
from ‘‘rigid’’ to fluid spheres. The region of ellipsoidal bubbles extends over 
different ranges of Reynolds numbers for the various liquids. The drag coef- 
ficients of spherical cap bubbles are independent of bubble size and have a 
constant value of 2.6. 

For bubbles (equivalent radius 0.03 to 0.30 cm) rising in tap water or 
in water containing certain surface-active substances, experiments show an 
increase in drag as compared to bubbles in pure water. 

Results of tests to determine the effect of the container walls on the 
velocity of rise are presented. A description of the experimental apparatus 
is given. A summary of the theoretical and experimental work of other in- 


vestigators is also included. 


INTRODUCTION 


The tests described in this report are a continuation of experiments given in a previous 
Taylor Model Basin report.! These experiments were initiated in conjunction with the work 
under projects NS 713-201 and NE 051-237. The present tests, continued under NS715-102, 
were conducted for the purpose of investigating the motion of bubbles rising under the 


lpeferences are listed on page 43. 


influence of gravity as an initial step for obtaining information on the behavior of bubbles in 
variable pressure fields. 
A body rising or falling under the influence of gravity reaches a constant velocity 


(terminal velocity) when all forces acting on it are in equilibrium: 
Drag + Buoyant force + Weight = 0 


For rigid bodies, the drag will, in general, be a complicated function of the geometry of the 
body, the velocity, and the physical properties of the medium, i.e., the density and viscosity. 
For fluid bodies, such as drops and gas bubbles, the function is further complicated by the 
fact that the body may be of changeable shape and that properties of the fluid inside the glob- 
ule, such as density and viscosity, and interfacial effects may also be important factors. In 
general, the shape that the fluid globule assumes is some complicated function of the hydro- 
dynamic, viscous, and interfacial forces. 

The drag of fluid bodies may either be equal to (as is the case for small bubbles) or 
less than that of the corresponding rigid body depending upon the conditions at the interface. 
In the former case there exists, effectively, a rigid surface at the interface; in the latter case, 
the fluid particles at the boundary have, in contrast to rigid bodies, nonvanishing tangential 
velocities. The circulation inside the fluid body thus reduces the drag of the body. 

The experiments described in this report consisted of the determination of the terminal 
velocity, shape, and path of single air bubbles rising freely in various liquids as a function 
of bubble size. The possible effect of the walls of the container on the velocity of rise of 
the bubble was also investigated. A summary of pertinent theoretical and experimental work 


of other investigators is included. 


THEORETICAL SOLUTIONS 


Theoretical solutions for the drag of rigid and fluid spheres, moving slowly in an infi- 
nite medium, have been obtained for the following boundary conditions at the surface of the 
sphere: 

1. Rigid spheres 
a. Stokes’ solution? 
(1) Velocity 
(a) The velocity vanishes. 
2. Fluid spheres 

a. Hadamard-Rybezynski’s solution?’3*4 

(1) Velocity 
(a) The normal velocity component vanishes. 


(b) The tangential velocity components at both sides of the surface are equal. 


(2) Stresses 
(a) The normal and tangential stresses at both sides of the surface are equal.* 
b. Boussinesq’s solution>’® 
(1) Velocity 
(a) The normal velocity component vanishes. 


(b) The tangential velocity components at both sides of the surface are equal. 


(2) Stresses 


(a) The normal stress at the inside of the surface is larger than the stress 
at the outside due to the dynamic surface tension.** 


(b) The tangential stress at the surface is increased across the surface due 
to the dynamic increment of surface tension. 


The drag of a sphere in an infinite medium of uniform velocity U thus becomes: 


1. Rigid sphere: D=6nrprtU 


2 3p’. 
2. Fluid sphere: a. D=6rprU Beal ase 2 
Bin Opn 


2 3’ 
b. D=67prU eshte pt 


e@+ 3r(pt pe’) 


where D is the drag, 
p is the coefficient of viscosity of the medium, 
r is the radius of the sphere, 
yz’ is the coefficient of viscosity of the fluid inside the sphere, and 


e is the coefficient of surface viscosity. 


Using the condition of equilibrium for a sphere rising under the influence of gravity, 


we obtain for the rigid case: 


4 selene, 
Cor Se) re pg 


2 9g 
U seg ran (p - p’) Stokes’ Law 
mn 


*The pressure increase across the surface due to surface tension (= 22) was neglected in Hadamard’s analysis. 
Inclusion of this pressure drop in the boundary condition for the normal stress does not change the results. That 
is to say, surface tension as manifested only in a pressure increase inside the fluid sphere does not affect its 
motion. (This result is also obtained by putting, in Boussinesq’s analysis, the coefficient of surface viscosity, 


see subsequent footnote, equal to zero.) 


**Boussinesq assumed nee S dynamic surface tension exists at interfaces in motion. Its magnitude is given by 
the sum of the usual (static) surface tension and the dynamic increment. The dynamic increment varies over the 
surface of the sphere and at a given point is proportional to the rate of dilatation at that point. The constant of 
proportionality is called (due to its similarity to the viscosity coefficient) the coefficient of surface viscosity. 


(Surface viscosity has the dimensions mass/time, while the dimensions for viscosity are mass/length x time.) 


where p is the density of the medium and p’ is the density of the sphere. 


For the fluid sphere: 


2 , 
u=2 TY - p’) 3y"+ 3p Hadamard-Rybczynski’s Law 
2 ip 3y’+ 2p 


2 , 
u=2 YY(o =p’) e+ 2r (utp) Boussinesq’s Law 
a 4 e+r(2y+ 3p’) 


For the case of bubbles, where 1.’ << » and p’ << p, the last two expressions reduce to? 


and 


In the last equation, the factor ee approaches 1 for r approaching zero or e very 
@+ 4Tu 


large; it approaches 3/2 for e << rp, i.e., for large r or for e approaching zero. Hence, 
for very small bubbles, Boussinesq’s solution approaches Stokes’ law as a limit, while the 
other limit is Hadamard-Rybezynski’s solution. 

From the boundary conditions as stated above, it is obvious that, for the Hadamard- 
Rybezynski and Boussinesq solutions, circulation exists inside the bubble. For the Stokes 


solution, of course, there is no circulation inside the sphere. 


DIMENSIONAL ANALYSIS 


Since an analytical solution for the drag of fluid bodies over a large range of sizes is 
hardly attainable, dimensional analysis of the phenomenon may serve to correlate the experi- 


mental results. 
In the case of fluid bodies the following physical variables are usually considered as 


pertinent: 


U Velocity of the body 

g Acceleration due to gravity 

p Density of the fluid medium 

p’ Density of the fluid body 

l Length parameter of the fluid body 


m Coefficient of dynamic viscosity of the fluid medium 


pe Coefficient of dynamic viscosity cf the fluid inside the fluid body 
oO Interfacial tension” 


The complete set of dimensionless products will therefore contain five such products. In 


principle, it is immaterial which complete set is chosen for the representation of the phenom- 


enon. For example, we may use 
in(Gp Re, We, ee 2) S 0 
wp 
or 


iis (Cp> Re, M, le", 12) =0 
p 


and so on, where Cp is the drag coefficient, 


Re is the Reynolds number i. zoe . 
pe 


2 
We is the Weber number (- 2H 4 


Oo 


Re’ is the Reynolds number inside the fluid body (- 2s.) and 


mn 


4 ay 
Mis a dimensionless parameter (-22 
3 
po 


If the density and viscosity of the gas inside a bubble are considered negligible, the physical 


variables are reduced to six. The dimensionless products then take the form of 


f,(Cp> Re, We) =0 


or 

f, (Cp, Re, M) =0 
or 

fs (Cp, We, M) =0 
etc. 


When experimental data on gas bubbles are plotted in terms of dimensionless products, 
complete correlation will be obtained provided the variables chosen are all the variables upon 
which the phenomenon depends. 

In the case of bubbles, it is most convenient to use a length parameter which is based 
on its volume rather than a physical dimension as is customary for rigid bodies. The length 


parameter chosen is the equivalent radius r, where 


_ 3/ volume 
: 2 


3, 


*The coefficient of surface viscosity is not included since there is no experimental evidence that dynamic 


surface tension, as postulated by Boussinesq, exists. 


For bubbles rising at their terminal velocity the drag coefficient C can then be written as 


(8/3) gr, 
D [tee -7- 
U2 


instead of its usual form of D/(1/2)p U2 A where A is the projected area. For the analytical 
solutions of Stokes and Hadamard-Rybezynski the drag coefficient becomes, respectively: 


== (Rigid sphere) 
16 : 
Cy =— (Fluid spher 
p=z2 (Fluid sphere) 


For the special case in which only four variables, namely the velocity, the accelera- 
tion of gravity, the density of the fluid medium, and the equivalent radius, are taken as per- 
tinent, only one dimensionless group, the drag coefficient, is obtained, i.e., 


Cp = constant 


This solution will be shown to apply to the region of very large (spherical cap) bubbles. 


PREVIOUS EXPERIMENTAL WORK 
RATE OF RISE OF BUBBLES 


Interest in the motion of air bubbles has existed for many years. The work on bubbles 
has, however, been mostly experimental in nature. Exceptions are an attempted theoretical 
analysis by Theremin’ in 1829, the analytical solution of Hadamard,?»3 Rybczynski,* and 
Boussinesq,°’® and the dimensional analysis of Schmidt® and Rosenberg. ! 

The early experimental work on bubbles was largely concerned with very small bubbles 
and was carried out for the purpose of determining the extent of Stokes region. Allen? deter- 
mined the rate of rise of air bubbles in water and in aniline up to bubble radii of 0.04 and 
0.06 cm, respectively. Arnold!° measured velocities of small air bubbles in olive oil and in 
aniline. Bond and Newton!! investigated air bubbles in syrup and in water glass (sodium 
silicate). 

The range of bubble sizes was extended by other investigators!?»!3 who were mainly 
interested in the problem in connection with air-lift pumps, !4-!7 gas absorption, 18-2° or 
propagation of sound in liquids.21 These experiments were carried out in water. In subse- 
quent years, some investigations were also made in liquids other than water. Davies and 
Taylor? used nitrobenzene as well as water and measured velocities of large bubbles. 
Temperley and Chambers?? extended the range of Taylor’s experiments in water to bubbles 
of equivalent radii up to approximately 6 cm. Bryn?* made tests in various water-glycerine 


and water-ethyl alcohol mixtures. Robinson?> measured the rate of rise of small air bubbles 


in lubricating oils. His results, however, show considerable scatter. Reports by Pickert,2® 
Pekeris,27 Worster,2® and Datta et al 2? give summaries of the results of experiments of 
other investigators. 

Gorodetskaya*° investigated the effect of surface-active substances on the rate of 
rise of air bubbles in water. Further tests on air bubbles are reported in References 31-36. 
In addition, a limited number of tests using gas bubbles of oxygen, nitrogen, and a mixture 
of carbon dioxide and oxygen have been carried out in artificial sea water.?7"38 A number of 
tests with oxygen bubbles were also conducted in water?’ and in aqueous solutions of sodium 
hydroxide.?9 Recently, Stuke?° investigated the rate of rise of oxygen bubbles in pure (pre- 
sumably distilled) water and in water containing surface-active substances. The results of 
the tests with gas bubbles (given in the Appendix) show no significant change in the rate 
of rise of the bubbles with change in the gas inside the bubble. 

Because of the scatter of previous results of experiments on the rate of rise of air 
bubbles in water, Rosenberg! repeated these tests for a large range of air bubble sizes. He 
showed the geometric similarity between large bubbles of spherical cap shape and suggested 
the use of three dimensionless parameters, the drag coefficient, the Reynolds number, and 


the parameter M for describing bubble motion in liquids. 


WALL EFFECT 


Previous investigations on the motion of air bubbles in liquids were, with a few ex- 
ceptions, conducted in containers of limited dimensions. Only for very large bubbles did 
Exner!3 and Bryn?* make their measurements in lakes. Inasmuch as the effect of the walls 
of the container on the rate of rise of bubbles was unknown, it was generally neglected. 
Miyagi !2 conducted a few tests in containers of different sizes and found that a reduction of 
4 percent in the rate of rise occurred for the range of bubble sizes investigated. Dubs*! de- 
rived, from energy considerations, an analytical expression for the wall effect and concluded 
that a bubble of the same radius as its cylindrical container has a velocity of rise of zero. 

It is clear that this conclusion is in error, as experiments on cylindrical bubbles have shown. 

As indicated previously, no analytical solution has yet been obtained for flows beyond 
the region of slow flow; consequently, the much more difficult problem of also including a 
finite boundary in the equations of motion becomes less capable of solution. For very slow 
flow about rigid spheres moving in an infinite cylindrical container, Ladenburg*? obtained 
an analytical solution for the effect of the boundary on the drag and, consequently, the veloc- 
ity of the sphere. McNown et al‘? arrived experimentally at a wall correction coefficient for 
rigid spheres descending in a cylindrical container. Since air bubbles of small volume rising 
in water behave essentially like rigid spheres, this correction factor may be applied to such 
bubbles as long as the flow is still in the Stokes region. 

With the exception of a number of tests 117 for large bubbles, no data concerning the 


effect of the boundary on the rate of rise of gas bubbles beyond the region of slow flow are 


available.’ However, Méller*4 showed, by means of a dye technique, that the flow about a 
rigid sphere at a lower Reynolds number and larger boundary dimensions was identical to that 
at a higher Reynolds number and smaller boundary dimensions and therefore that the effect 

of the walls was to stabilize the flow about the sphere. These results for rigid spheres be- 
yond the Stokes region of flow at least suggested the possibility of a similar effect for the 


motion of gas bubbles. 


CYLINDRICAL BUBBLES 


At this point, it may be of interest to mention experiments on a special form of finite 
boundary dimensions, i.e., the case of cylindrical bubbles. This term was first used by 
Gibson*® and applies to the type of bubble formed when a long cylindrical tube filled with 
liquid is emptied from below or when a large amount of air is introduced through the bottom 
of the tube. 

Gibson investigated the velocity and shape of these bubbles in water. Ward and Kess- 
ler!® conducted tests in pipes of various diameters. Hattori*® was interested in the problem 
in connection with the possibility of evaluating the surface tension of a liquid. Hence, he 
was concerned with tubes of smal! diameter, since the so-called critical tube diameter (below 
which the bubble no longer rises but remains stationary) is a function of the surface tension 
of the liquid. Dumitrescu4” obtained an analytical expression for the velocity of a cylindrical 
bubble by neglecting viscous and surface tension forces, thus reducing the problem to one of 
potential flow. The differential equation for the velocity potential together with the existing 
boundary condition yields a solution for the velocity of rise as a function of the tube diameter 
only. His experimental tests in water show that for a tube of sufficiently large diameter 
(3 cm for water at room temperature) the measured velocities agree very closely with the 
theoretical values. Therefore, for large bubbles, the physical properties of the liquid no 
longer have any effect on the flow about the bubble and the bubbles are geometrically similar. 
Davies and Taylor2? investigated the shape and rate of rise of cylindrical bubbles in order 
to obtain a better understanding of the pressure distribution of spherical cap bubbles in an 


infinite medium. 


’ SCOPE OF THE PRESENT INVESTIGATION 


The present investigation was initiated in connection with a program of study of the 
behavior of air bubbles in water at variable pressure gradients. Since extensive experimenta- 


tion was required for direct experimental study of the motion of air bubbles in such pressure 


*A pap ~ by Coppock and Meiklejohn>4 has recently come to the attention of the authors. From tests con- 
ducted with air bubbles in water, they conclude that no wall effect exists for bubbles ranging in equivalent 


radius from 0,01 to 0.1 cm rising in a tube of 5 cm diameter. 


fields,” an alternate. approach appeared more feasible. It consists of calculating the motion 
of the bubbles in variable pressure fields from a knowledge of bubble drag at various constant 
pressure gradients. Experimental data on the drag of air bubbles at various pressure gradients 
are essential in this procedure. However, only data on the motion of bubbles in pressure gra- 
dients produced by gravity were available. Therefore, information on bubble motion in water 
at pressure gradients other than gravity became necessary. This information could be obtained 
by investigating the rise of bubbles in liquids having the same physical properties as water 
with the exception of the density. With all other properties of the liquids identical, varying 
the density would be equivalent to varying the pressure gradient. This approach, however, is 
not practicable, since there are no liquids available which possess such properties. The 
other approach is to investigate the rise of air bubbles in various liquids having different 
physical properties and then to attempt to correlate the results in terms of nondimensional 
parameters. The available information on the rise of bubbles in different liquids was too 
meager to allow definite conclusions regarding the significance of the parameters suggested 
by Reference 1. The present investigation was therefore initiated with the purpose of deter- 
mining the nondimensional parameters for bubble rise by investigating bubble motion in a 
number of liquids of different physical properties. If it were found that the motion of air bub- 
bles rising freely in a liquid, that is to say the motion in the pressure field produced by gravi- 
ty, could be described, for example, in terms of the drag coefficient, the Reynolds number, 
_and the parameter M(= gu*/po*),** the results thus obtained could be used in evaluating the 
drag of bubbles in water at pressure gradients other than that produced by gravity. To do 

this it would have to be shown that the nondimensional parameters used for the freely rising 
bubbles are also applicable to other pressure fields. This might be accomplished, for example, 
by comparing the results of a bubble experiment in water in a nongravity pressure gradient 
with the results of bubbles rising freely in various liquids at identical ‘‘¥’’? number. 

By conducting the tests on the rise of bubbles in various liquids in a large tank, the 
possibility of the effect of the tank walls on the velocity of the bubbles is eliminated. Since 
the high cost of many desirable liquids makes the use of a large tank impractical, it became 
necessary to determine the possible effect of the walls on the velocity over the range of bub- 


ble sizes to be tested. This investigation consequently acquired two purposes: 


1. The determination of the effect of variation of liquid properties on the motion of air 
bubbles. 


2. The evaluation of wall effect. 


*Exploratory experiments of such a nature are reported in Reference 35. 


**This parameter is given in a more general form as (ut Vp)/(p20°) (where Vp is the pressure gradient; for a 
gravity field Vp = p&). Therefore, for a specific liquid, it is proportional to the pressure gradient. 
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EXPERIMENTATION 


The experimental study consisted of measuring the terminal velocity of individual bub- 
bles of various sizes rising in eight liquids. It also included the determination of the effect 
of the walls of the container on the bubble velocity. The experimental apparatus consisted 
of test tanks, means for determining the physical properties of the liquids, a regulated bubble 
supply, and means for measuring bubble size and velocity. Details of the experimental appa- 
ratus and procedure, the generation of the bubbles, and the test liquids are given in the fol- 


lowing paragraphs. 


TEST TANKS 


In order to obtain free bubble rise and to reduce the effects of such boundaries as the 
bottom of the tank and the free liquid surface, the containers for the liquids had to be at 
least 2 ft high. The tests were performed in three transparent wall tanks; the large one was 
of 3 x 3 ft cross section and 5 ft height, the medium one was of 1 x 1 ft cross section and 
3 1/2 ft height, and the small one was of 6 x 6 in. cross section and 2 ft height. In addition, 
tests were also performed in an insert of 6 x 6 in. cross section and 20 in. height, placed in 
the center of the medium tank. 

Since the large tank was of sufficiently large dimensions, no significant wall effect 
was expected. The dimensions of the medium tank were chosen large enough to reduce wall 
effects, yet small enough to allow use of a variety of liquids. The small tank and the insert 
provided an additional tank size. It was intended that if wall effect existed, the results ob= 
tained in the finite containers would be extrapolated to the case of an infinite medium. 


TEST LIQUIDS 


The eight test liquids used were water (at three different temperatures), Varsol, 
methyl alcohol, turpentine, water containing 0.42 percent (by volume) Glimsaes mineral oil, 
and two corn syrup-water mixtures. Turpentine was selected as one of the test liquids be- 
Cause at room temperature it has the same viscosity as cold water. One of the corn syrup 
mixtures had approximately the same viscosity as the mineral oil (see Table 1). 


The viscosity of the liquids was measured by means of ordinary and modified Ostwald 
viscosimeters. The accuracy of measurement of viscosity was 1.5 percent and 0.5 percent, 


*A few of the previous Taylor Model Basin tests 1 


were repeated in the large tank to observe any change in 
results. These previous tests were conducted in a tank of 4 1/2 X 25 ft cross section and 9 ft height with 8 ft 


depth of filtered water at room temperature, using one end of the tank for the tests. 
**A trade name (Standard Oil Company) for mineral spirits (heavy naphtha), a petroleum distillate. 


***Glim (Antarox A-480), a surface-active agent, is the trade name (B.T. Babbitt, Inc.) of a non-ionic, liquid 
detergent, a condensation product of ethylene oxide and lauryl alcohol. 
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TABLE 1 
Summary of Liquid Properties 


Surface 
Temperature Viscosity Density Tension 
Liquid im p o 
deg C poises gm/cc dynes/cm 


Water 0.26 x 10°20 
Water 0.24 x 19710 
Cold Water HN) 52 yee 
Hot Water 0.307 x 1071! 
Glim Solution 2.18 x 102° 
Mineral Oil 1.45 x 10-2 
Varsol 4.3 x 10-10 
Turpentine 24.1 x 19710 
Methyl Alcohol 0.89 x 19710 


62 percent Corn Syrup 0.155 x 10°3 
and Water : 


68 percent Corn Syrup 0.212 10-2 
and Water $ 


56 percent Glycerine ! 57 See] One 
and Water (Bryn2*) 


42 percent Glycerine , 4.18 x 1078 
and Water (Bryn2*) 


13 percent Ethyl Alcohol 117 x 10° 
and Water (Bryn2*) : 


Olive Oil (Arnold!) 0.716 x 1072 
Syrup (Bond?!) 0.92 x 10° 


respectively. The surface tension was determined by the capillary-rise method (accuracy of 
measurement: 3 percent) and the specific gravity of the liquids was obtained by means of 
hydrometers (accuracy of measurement: 1 percent). These physical properties were measured 
following the completion of each test. They are summarized in Table 1 together with those of 
liquids used by several other investigators. 19 11,24 

For Varsol and water (room temperature and hot), tests were conducted in all three 
tanks; for cold water and mineral oil, in the medium tank and insert; and for all other liquids, 
in the small tank only. 
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Camera 


i) 
Figure 1 - Experimental Setup 


BUBBLE GENERATION 


Small bubbles were generated by means of hypodermic needles and glass nozzles of 
various sizes. The larger bubbles were obtained by use of a dumping cup, which was inverted 
to release the air bubble. The nozzles and needles were connected to a brass tube which 
was fastened to a sliding mechanism (Figure 1). This sliding mechanism allowed the tips of 
the various nozzles to be placed at the identical position. The air was supplied from a com- 
pressed air bottle. A needle valve regulated the air flow so that bubbles were released at the 


interval desired. 


DETERMINATION OF BUBBLE SIZE 


The bubble size was determined by ‘‘weighing’’ a sufficient number of bubbles in the 
inverted funnel (Figure 1) by means of an analytical balance. Since the density of air is 


negligible in comparison to that of a liquid, the difference in balance reading equals the 


*This device eliminated focusing of the camera after each change of nozzle. The camera was used to deter- 


mine velocity, path, and shape of the bubble. 
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buoyancy of the bubbles (i.e., it equals the volume of the bubbles times the density of the 
liquid). The change in balance reading was always at least 0.2 gm, resulting in an accuracy 
of measurement of volume of 1 percent. The volume of the individual bubble was obtained by 
dividing the total volume by the number of bubbles collected in the funnel. A comparison of 
photographs of different bubbles showed that the bubble size did not vary if the frequency of 
bubble generation remained constant. Large bubbles from the dumping cup were weighed 
individually. The volume of the bubble was adjusted for the change in pressure due to differ- 
ence in depth between the level at which the rate of rise is determined and the level of the 
inverted funnel. This was done by use of the general gas law at constant temperature, taking 
into account the partial pressure of the saturated vapor at test temperature. - 

Tiny spherical bubbles could not be generated at a frequency to allow a sufficient 
number to be collected in the funnel, hence their size was determined from the photographic 


record. No correction for change in depth is then needed. 


OTHER EXPERIMENTAL TECHNIQUES 


To avoid any changes in the volume of the bubble due to air interchange with the 
liquid, the latter was saturated with air prior to actual testing. This was accomplished by 
stirring the liquid and by blowing air through it. 

The temperature of the liquids, with the exception of cold and hot water, was room 
temperature, which varied little throughout the day. Water was cooled by circulation through 
a water cooler; it was heated with immersion heaters or obtained directly from the hot water 
faucet. Both filtered and tap water were used in the tests. 

Uniformity of liquid temperature was achieved by means of mechanical stirring before 
each test. Frequent checks of temperature at various locations inside the tank were made by 
means of immersion thermometers. In the process of stirring, small bubbles appeared in the 
liquid. The irregularity of motion of these small bubbles, which were still present after com- 
pletion of the stirring, served as an indication of the presence of residual turbulence in the 
liquid. In sufficient time, the motion of the small bubbles always became regular and hence 


indicated that the residual turbulence, if still present, was not large enough to affect the 


*Details of this correction are as follows: a 


1 0) 


Py-P 

where Vy is the volume of the bubble at the camera level, 
Vo is the volume of the bubble as determined by weighing, 
Po is the absolute pressure at the funnel, 


Py is the absolute pressure at the camera level, and 


Pp. is the vapor pressure of the liquid at the test temperature. 
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motion of the bubbles. The actual tests were not begun until all the small bubbles reached 
the surface of the liquid. 

The rate of bubble flow was then regulated by the needle valve so as to release bubbles 
with a minimum spacing of 24 in. This reduced the effect of the wake created by the passage 
of a bubble on the motion of a bubble following.” At higher bubble rates, the velocity of the 
individual bubble is increased. 

The same precaution was observed for the larger bubbles that were formed by dumping. 
An additional precaution was to rotate the dumping cup with steady speed in order to avoid 
splitting of the bubble or the formation of satellites upon release. The slow passage of the 
air through the brass tube inside the tank allowed the air to reach the temperature of the li- 
quid. Contact of the air at the nozzle tip or inside the dumping cup with the liquid allowed 
saturation of the air with liquid vapor, so that the air bubble can, in each instance, be as- 


sumed to be saturated with the vapor of the liquid in which it rises. 


MOTION PICTURES AND THEIR EVALUATION 


The velocity, path, and shape of the bubbles were obtained from motion pictures made 
with a Mitchell 35mm camera using a special lens attachment to permit close-ups. Film speeds 
of 25 to 35 frames per second and back lighting from a white reflector were used. For the first 
few tests, the film speed was obtained by photographing a rotating clock dial; subsequently a 
neon timing light with a 60-cycle voltage source was utilized. The film speed was determined 
from the marks of the timing light on the film. The field of the camera varied from 1.4 x 1.8 
to 1.75 x 2.3 in. depending upon the refractive index and horizontal depth of liquid. A trans- 
parent scale photographed in the plane of the bubble provided the distance scale factor for 
the evaluation of displacement and size. The camera lens was placed at approximately the 
midpoint between the liquid level and the bottom of the tank for all tests. The camera loca- 
tion was in each instance sufficiently above the nozzle tip so that the bubbles reached their 
terminal velocity before passing in front of the camera. A summary of camera location and 
depth of liquids in the tanks is given in Figure 2. 

Changes in bubble volume due to differences in liquid depth were minimized by making 
velocity measurements over a very short vertical displacement (less than 2 1/2 in.). The 
rate of rise of bubbles was determined by measuring the displacement of a bubble from a 
reference point on successive frames of the film by means of a Bausch and Lomb contour- 
measuring projector using a magnification of twenty-five (Figure 3). These displacements 
were then plotted against the frame number. The straight-line plot indicates that the veloc- 
ity of the bubble remained constant during the time it passed the field of the camera. From 
the slope of the line, the frame speed, and the scale factor, the velocity of the bubble is 


computed. 
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*Napier~- showed the absence of proximity effect for air bubbles in water, ranging in equivalent radius from 


0. 14 to 0.38 cm, if the frequency was below 30 bubbles per minute. 
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RESULTS 
TERMINAL VELOCITY OF AIR BUBBLES 


The results of tests to determine the velocity of rise of air bubbles in various liquids 
are most cenveniently presented as a function of the equivalent radius of the bubble, defined 
as the radius of a sphere having the same volume as the bubble. Figures 4-13 show the ter- 
minal velocity of air bubbles rising freely in tap (unfiltered) and in filtered water (including 
data from other investigators), in water containing Glim, in mineral oil, Varsol, turpentine, 
methyl alcohol, and two corn syrup-water mixtures as a function of the equivalent radius. 
Figure 14 presents Bryn’s results in an ethyl alcohol-water mixture and two glycerine-water 
mixtures. Figure 15 summarizes all velocity curves (except those for tap water). A compila- 
tion of the properties of the liquids is given in Table 1 (see page 11). 

In general, the results as seen from Figure 15 indicate that for small (spherical) air 
bubbles of given volume, the viscosity of the liquid is the most important property determining 
the rate of rise. Very large bubbles (spherical caps) rise independently of the properties of 
the liquid. 


WALL AND PROXIMITY EFFECTS 


As indicated previously, the effect of the container walls on the velocity of a bubble 
had to be determined if the results of tests conducted in a tank of limited dimensions were 
to be applied to bubble motion in an infinite medium. Tests were, therefore, conducted in 
tanks of different sizes in water, Varsol, and mineral oil. Figure 4 gives the results of tests 
conducted with filtered water by several investigators including the Taylor Model Basin. 
The cross sections of the containers used are also indicated in the figure. No wall effect 
is noticeable from these results. For example, Gorodetskaya’s results for bubbles ranging 
from 0.01 to 0.07 cm rising in a tube of 5 cm diameter show no wall effect when compared 
with results of tests conducted in larger containers. The results of the present experiments, 
given in Figures 5, 6, 7, 9, and 10, show within experimental accuracy, the absence of any 
wall effect for the range of bubble sizes tested. Subsequent tests in the other liquids were 
made in the small tank only and the results may be applied to the case of an infinite medium. 

No systematic investigation was made of vertical proximity effect, i.e., the effect of 


the wake created by the passage of a bubble on the motion of a bubble rising at a distance 


*The results of the 81 percent (by weight) glycerine-water mixture have been omitted. Bryn presented these 
results in terms of drag coefficient and Reynolds number, from which the terminal velocity can be computed. 
In the region of bubble size where the rate of rise is shown to be a function of size only (hence a common veloc- 
ity curve for all liquids; see Figure 15), the velocity curve for the 81 percent mixture falls appreciably above 
the common curve. The discrepancy is probably due to erroneous evaluation of the two dimensionless parameters 


for the 81 percent glycerine-water mixture. 
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Figure 4 - Terminal Velocity of Air Bubbles in Filtered or Distilled 
Water as a Function of Bubble Size 
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Figure 5 - Terminal Velocity of Air Bubbles in Tap Water 
as a Function of Bubble Size 
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Figure 6 - Terminal Velocity of Air Bubbles in Cold Filtered Water 
as a Function of Bubble Size 
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Figure 7 - Terminal Velocity of Air Bubbles in Hot Tap Water 
as a Function of Bubble Size 
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Figure 8 - Terminal Velocity of Air Bubbles in Water Containing Glim 
as a Function of Bubble Size 
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Figure 9 - Terminal Velocity of Air Bubbles in Mineral Oil 
as a Function of Bubble Size 
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‘Figure 10 - Terminal Velocity of Air Bubbles in Varsol 
as a Function of Bubble Size 
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Figure 11 - Terminal Velocity of Air Bubbles in Turpentine 
as a Function of Bubble Size 
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Figure 12 - Terminal Velocity of Air Bubbles in Methyl Alcohol 
as a Function of Bubble Size 
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Figure 13 - Terminal Velocity of Air Bubblés in Corn Syrup-Water Mixtures 
as a Function of Bubble Size 
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Figure 14 - Terminal Velocity of Air Bubbles in Liquids as Obtained 
from Bryn’s Data 
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Figure 15 - Comparison of Terminal Velocities of Air Bubbles in Various Liquids 
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below. Such effects were avoided in the experiments by sufficient spacing between the bub- 
bles. However, the results of a few special observations indicate that proximity effects may 
be appreciable. For example, tests in mineral oil show an increase of 9 percent and 39 per- 
cent for bubbles of equivalent radius of 0.17 cm, rising 7.7 cm and 3.2 cm apart, respectively. 
Napier?° observed an increase of 6 percent for bubbles of 0.14 cm radius in water, rising 

6 cm apart. The presence of the wake in the liquid thus results in higher velocities of rise 
of the bubble. 


NONDIMENSIONAL PRESENTATION OF BUBBLE DATA 


In a previous section it was pointed out that presentation of the experimental data on 
air bubbles in terms of dimensionless products gives complete correlation provided the vari- 
ables considered in the analysis are complete and pertinent. The results of the Taylor Model 
Basin bubble tests and those of Arnold,!° Bond and Newton,!! and Bryn?* are given in terms 
of the drag coefficient, Reynolds number and the parameter M in Figure 16.* Figure 17 pre- 
sents the bubble data in terms of the drag coefficient, Weber number, and the parameter M. 
The curve for filtered or distilled water at a temperature of 19 deg C was drawn through points 
obtained from the experiments of Bryn and the Taylor Model Basin tests. 

Examination of Figure 16 or 17 shows no systematic arrangement of the curves with 
change in the parameter M, which is constant for a specific liquid. It can, therefore, be con- 
cluded that neither of the nondimensional sets presented nor any other complete set using 
the same six variables (namely, velocity, acceleration of gravity, density and viscosity of 
the liquid, surface tension, and equivalent radius) is sufficient for a complete description of 
bubble motion. : 

The question now arises whether correlation of bubble data could be obtained by using 
two additional dimensionless parameters, for example, the liquid to air viscosity and density 
ratios or the Reynolds number inside the bubble and the density ratio, etc. The results of 
the experiments conducted do not permit conclusions regarding the importance of these 
parameters. A short discussion of the significance of the internal Reynolds number will be 


given in subsequent sections. 


SPHERICAL BUBBLES 


It was observed in the experiments that, as the bubble size was increased, a change 
in bubble shape from spherical to ellipsoidal to spherical cap shape occurred in all liquids. 
Very small bubbles are spherical. Larger bubbles are flattened, i.e., ellipsoidal in shape, 


whereas very large bubbles assume a spherical cap shape. Of course, the volumes at which 


*The results for tap water and for water containing Glim are not shown. They will be discussed in subsequent 


sections. 
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these transitions occur vary with the liquids. Photographs of typical shapes are shown in 
Figure 18. It should be noted here that some of the shapes shown in these photographs are 
instantaneous shapes, since the shape of large bubbles does not remain constant during the 
ascent. An exception are bubbles rising in a highly viscous medium (e.g., mineral oil and 
corm syrup). 

The results for the spherical bubbles only are plotted in terms cf the drag coefficient 
and the Reynolds number in Figure 19, with Reynolds numbers ranging up to about 400. The 
drag curve for rigid spheres is also included.*® From it, the following can be observed: The 
drag curves of spherical bubbles in the various liquids fall between two limiting curves. As 
upper limit, the drag curve of rigid spheres is obtained, while the lower limit is the drag curve 
for fluid spheres. With decreasing Reynolds number, the rigid sphere curve connects with the 
straight line of Stokes’ Law, while the fluid sphere curve connects with the line of Hadamard- 
Rybezynski’s Law. The curve for the fluid spheres was obtained by drawing the lower en- 
velope to the experimental curve; its accuracy can be confirmed by additional tests in other 
liquids or by extension of the theoretical solution into regions beyond that of very slow flow. 

It will be noted from Figure 19 that the curve for mineral oil, for example, follows the 
straight line of Hadamard-Rybezynski’s Law over a certain region of Reynolds numbers. 

This indicates that the boundary conditions assumed in the analytical solution for fluid 
spheres are actually fulfilled and that circulation exists inside the bubble. Circulation in- 
side bubbles has been observed experimentally.*9 

The experimental curves of Figure 19 also indicate an interesting aspect of the phe- 
nomenon of bubble motion, namely that with decreasing Reynolds number, the drag coefficient 
of the bubbles becomes equal to the drag of rigid spheres. This transition may occur at a 
Reynolds number of about 40 (as for filtered and distilled water) or may not take place until 
very low Reynolds numbers are reached, i.e., well within the region of slow flow (as for olive 
oil!® or very viscous syrup!?). Thus, from the experimental data available, it appears certain 
that tiny air bubbles rising in any liquid follow Stokes’ Law. 

For bubbles behaving like rigid bodies, thus indicating absence of motion inside the 
bubble, the internal Reynolds number (although nonvanishing) is of no significance in de- 
scribing the rising motion of the bubbles. Likewise, the internal Reynolds number cannot be 
used to predict the transition point at which the drag of the bubbles becomes less than that 
of corresponding rigid spheres. Beyond this transition point, the internal Reynolds number 
might be of importance in describing the motion of the bubbles. 

Surface tension tends to make the’surface area of the bubble as small as possible. 

For a given volume, the configuration of minimum surface area is a sphere. This effect of 


surface tension would be most pronounced for bubbles of small radii. 
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ELLIPSOIDAL BUBBLES 


‘For larger bubble sizes, the surface forces, which are essential in maintaining the 
spherical shape of a bubble, become smaller in comparison to the viscous and hydrodynamic 
forces, and flattening of the bubble occurs. This flattening to approximately an oblate 
spheroid results in higher drag as compared to a sphere of the same volume. Figure 20 shows 
the drag curves of ellipsoidal and spherical cap bubbles in terms of the Reynolds number. 
The estimated extent of the regions of ellipsoidal and spherical cap bubbles is indicated 


in the figure. It will be observed that the region of ellipsoidal bubbles for the various liquids 
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Figure 18 - Typical Shapes of Bubbles of Several Volumes in the Various Liquids 


occurs at different ranges of Reynolds number, that for liquids of low ‘‘M@’? number a minimum 
in the drag curves is reached at Reynolds numbers of the order of 250, and that these minima 
occur near the transition from spherical to ellipsoidal shape. Such minima are not obtained 
for liquids of high ‘‘M’’ number. The drag coefficients of bubbles in such liquids decrease 
until a constant value for the drag coefficient (spherical caps) is attained. In the ellipsoidal 
region the curves are arranged according to the magnitude of the ‘‘y’? number, indicating that 


the liquid properties contained in this parameter, namely surface tension, viscosity, and 
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density, are of primary importance in the motion of these bubbles. 

Thus, in this region the Reynolds number inside the bubble is of no importance in the 
description of bubble rise, since correlation was obtained in terms of drag coefficient, Rey- 
nolds number, and ‘‘’’ number. 

In Figure 21 the drag coefficient of the bubbles was plotted as a function of Weber 
number, [tis seen that, for the liquids tested, transition to a constant value of drag coeffi- 


cient (spherical cap region) is reached at a Weber number of about 20. 


SPHERICAL CAP BUBBLES 


When the viscous and surface tension forces become small relative to the hydrodynamic 
forces, the shape assumed by the bubbles is that of the so-called spherical caps. Typical 
shapes of these bubbles are shown in Figure 18. The upper surface is essentially spherical, 
while the lower surface varies from a highly irregular one for liquids of low viscosity to a 
smooth surface for very viscous liquids. The configuration of the upper surface results al- 
most exclusively from the hydrodynamic forces.” 

The geometric similarity of these bubbles was shown by Rosenberg! and Davies and 
Taylor,22 who determined a constant drag coefficient of 2.6 for them. The results of the 
present tests in a number of liquids confirm this value; Figure 21. The velocity of spherical 
cap bubbles of given size rising in any liquid can be determined from the constant value of 
the drag coefficient or directly from the velocity curve (Figure 15). For Cp = (8/3) gt -/U 


= 2.6, we obtain for the rate of rise of the spherical caps in all liquids 


U=1.02Ygr, 


Thus, the velocity of rise of these bubbles is a function of the bubble size only and 


not of the properties of the liquid (see special case of ‘‘Dimensional Analysis’’). 


PATH OF BUBBLES 


Figures 22-25 show representative paths and corresponding shapes of bubbles in the 
various liquids. Three types of motion of the bubbles were observed in the experiments: 
(1) rectilinear motion, (2) motion in a helical path, and (3) rectilinear motion with rocking. 
The motion of spherical bubbles is either rectilinear or helical. For ellipsoidal and spherical 
cap bubbles, all three types of motion can occur. It appears that the type of motion may be 
_ predicted from the value of the Reynolds number at which the motion takes place. Below 
Reynolds numbers of about 300 the motion is rectilinear. With increase in Reynolds number 
spiraling begins and increases in amplitude and frequency until a maximum is reached. At 


Reynolds numbers of about 3000, the spiraling disappears and only rectilinear motion with 


*See, e.g., Reference 22 
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rocking is obtained. For the bubbles rising in mineral oil and the corn syrup mixtures, only 
rectilinear motion without rocking was observed. The maximum Reynolds number reached 
during those tests was 150. 

The helical path of the bubbles assumes either a clockwise or counterclockwise direc- 
tion, depending upon conditions at generation. The velocity of rise of these bubbles is not 
affected by the sense in which the bubble revolves. The major axis of ellipsoidal bubbles 
is always perpendicular to the direction of motion. 

The oscillatory motion of bubbles is probably caused by the periodic shedding of vor- 
tices behind the bubble. Such vortex shedding has been observed experimentally for rigid 


spheres at the same magnitudes of Reynolds numbers as for bubbles.5° 


BUBBLES IN FILTERED AND TAP WATER 


The results of the experiments in filtered water at room temperature and incold filtered 
water are presented in terms of drag coefficient and Reynolds number in Figure 26. The re- 
sults of Bryn,?* Allen,? and Gorodetskaya?® are also included. It is seen that the drag curves 
at the two different temperatures coincide in the spherical and spherical cap region. In the 
region of ellipsoidal bubbles, the drag coefficient at a given Reynolds number increases with 
increase in ‘‘M’’ number. Gorodetskaya and Allen conducted their experiments in distilled 
water. A comparison of their experimental data with those obtained in filtered water shows, 
within experimental accuracy, no difference in the drag of air bubbles rising in filtered and 
distilled water. 

The drag coefficients for air bubbles rising in tap water at two different temperatures 
are given in Figure 27. Gorodetskaya’s results in tap water at room temperature are included 
in the figure. For comparison, the drag curves for bubbles in filtered water at room tempera- 
ture and for rigid spheres are also shown. Again, in the region of spherical and spherical cap 
bubbles, the drag curves at the two temperatures coincide. The value of the minimum drag 
coefficient is, however, greater than that of the corresponding filtered water. In general, for 
Reynolds numbers up to about 300, the drag curves of bubbles in tap water follow closely the 
curve of rigid spheres. 

Thus, the results of the experiments (given in Figures 5, 26, and 27) show that for bub- 
bles (ranging in equivalent radius from 0.035 to 0.25 cm) it is important whether the motion 
occurs in filtered (distilled) or tap water. 

In view of the fact that merely filtering the water was sufficient to produce a change in 
the drag of the bubbles, it is indicated that the presence of minute particles causes this 
change. Minute particles, most of which can be removed by filtering, are known to exist in 


ordinary tap water. Specifically, if such particles are present in the water a high concentration 


*It should be noted here that the physical properties of tap water did not differ from those of filtered water. 
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of these particles would be found at the surface of the bubble. The particles at the surface 
would travel with the bubble, hence imparting, in effect, a rigid surface to the bubble. With 
increasing bubble velocities, the shear forces become large in comparison to the forces hold- 
ing the particles to the surface and hence at a certain critical velocity no high concentration 
of particles on the surface can exist. 

Figure 28 compares the drag curve of bubbles in tap water obtained in the Model Basin 
study and the experimental data of air bubbles rising in water at room temperature as obtained 
by various other investigators. Reports of these other studies give no information regarding 
the type (tap, filtered, or distilled) of water used, but presumably these experiments were 


also conducted in tap water. 


EFFECT OF SURFACE-ACTIVE SUBSTANCES 


The effect of surface-active substances on the rate of rise of air bubbles has previous- 
ly been investigated by Gorodetskaya,°° who added small concentrations of various alcohols 
to water and concluded that, beyond a certain critical concentration of the surface-active sub- 
stance, the rate of rise of the air bubbles is not affected. Stuke?° ran experiments with oxy- 
gen bubbles rising in water containing small concentrations of caproic acid. The concentra- 
tions of the alcohols and the caproic acid were relatively small, hence the decrease in the 
surface tension was only about 1 dyne/cm. In the present study, the authors conducted tests 
in water containing Glim, a liquid detergent. The concentration of Glim (0.42 percent by vol- 
ume) was high enough to decrease the surface tension by 40 dynes/cm. No measurable change 
in the viscosity and density of the test liquid due to the presence of Glim was noted (see 
Table 1). This was also true for the alcohol and caproic acid solutions. Results from these 
experiments are presented in terms of the drag coefficient and Reynolds number in Figure 99." 

The drag curve for bubbles in the Glim solution, as well as the experimental data from 
the other investigators in water containing at least the critical concentration of the surface- 
active substance, follows the drag curve of rigid spheres to a Reynolds number of about 9200.* * 
In the region of Reynolds numbers of 10 to 200, the drag curve for bubbles rising in a pure 
liquid having an ‘‘M’? number very close to that of the Glim solution” follows the drag 
curve of fluid spheres. Thus, the motion of bubbles in water containing surface-active ma- 
terials cannot be compared with that of bubbles in pure liquids on the basis of drag coefficient, 


Reynolds number, and ‘‘’’ number, even in the region of ellipsoidal bubbles. Although the 


*As shown in the Appendix, no significant difference in the terminal velocity of oxygen and air bubbles rising 
in distilled water is obtained. Hence, inclusion of the results of the tests with oxygen bubbles in water contain- 


ing surface-active substances is justified. 


**For concentrations below the critical, the drag curve lies between that of pure water and the curve shown in 
Figure 29.40 


***The ‘?? number of Glim was 2,78 X 10° 10° that of Varsol, for example; was 4,3 X 10719, 
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range of bubble sizes did not extend fully into the region of spherical cap bubbles, it is quite 
certain that the presence of surface-active substances will not alter their rate of rise, which 
was shown to be independent of all physical properties of the liquid. 

The difference in behavior of the air bubbles must be sought in the behavior at the 
surface. A high concentration of molecules of the surface-active substances will be found 
at the surface of the bubble. As for the case of the tap water, these molecules would travel 
with the bubble and would impart, in effect, a ‘‘rigid’’ surface to the bubble, that is to say, 
would impose the condition of zero velocity at the boundary. 

Thus, the results of the tests given in Figure 25 show that these surface-active sub- 
stances increase the drag of the bubbles (in the region of bubbles having equivalent radii of 
0.03 to 0.30 cm); beyond the critical concentrations, any increase in concentration has rela- 


tively little influence on the drag of the bubbles. 


BUBBLES AS RIGID BODIES 


In previous sections it was shown that bubbles rising in pure liquids” behave essen- 
tially like fluid bodies over a large range of bubble size, but that below a certain critical 
size (the size being different for various liquids), the bubbles behave like rigid bodies, that 
is to say, the drag of the bubbles equais that of corresponding rigid bodies. 

A possible explanation of the anomaly of behavior of the gas bubbles is given by 
Boussinesq’s dynamic surface tension.> As pointed out under ‘‘Theoretical Solutions,” 
for small bubbles the effect of the dynamic increment of surface tension increases the drag 
to the value of corresponding rigid bodies. With increase in bubble size, this effect becomes 
negligible and the drag of the bubble equals that of a fluid body. The dynamic increment in- 
cludes a constant of proportionality (surface viscosity) which is a function of the two fluids 
composing the interface. Therefore, for air bubbles, the transition region from ‘‘rigid’’ to 
fluid bodies would be different for the various liquids. There is, however, no experimental 
evidence that dynamic surface tension, as postulated by Boussinesq, exists. 

From a ‘‘hydrodynamic’’ point of view, the reason for the transition of the bubbles to 
“rigid”? bodies is not clear. As mentioned earlier, the mere inclusion of surface tension as 
a pressure drop in the boundary conditions does not alter the analytical solution for fluid 
spheres. Hence, it appears that the presence of surface tension should have no effect on the 
motion of the bubble, except in maintaining the spherical shape. Thus the anomalous behavior 
of the bubbles cannot be explained in terms of ‘‘hydrodynamics,’’ but must be sought in terms 
of a surface phenomenon. If it could be shown analytically (if only in the region of slow flow) 
that equality of drag of corresponding rigid bodies and bubbles alsc implies equality of bound- 


ary conditions at the surface, then, as in the case of rigid bodies, the velocity of the entire 


*Mixtures such as the 13 percent ethyl alcohol-water mixture are included inthis category. 
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surface of the bubble must vanish. The surface must then be able to hold molecules of the 


pure liquid, just as in the case of tap water and surface-active substances the surface at- 
tracts and holds a high concentration of particles or molecules of the surface-active sub- 


stances. The molecules at the surface would travel with the bubble and hence would, in ef- 
fect, give the same boundary conditions as a rigid surface. As the shear forces become 
larger in comparison to the forces holding the molecules at the surface, ‘‘rigidity’’ at the 
surface cannot be maintained; circulation inside the bubble ensues and the drag of the bubble 
becomes smaller as compared to that of a rigid body. 

In addition, since it was not possible to correlate the results of the experiments on 
the motion of bubbles in the gravity pressure field in terms of nondimensional parameters 
formed from the usual liquid properties (viscosity, surface tension, density), further work on 
freely rising bubbles is necessary before the results obtained from such tests can be utilized 
and before the more complicated behavior of bubbles in variable pressure gradients can be 
understood. Particularly, an understanding of the reason for the transition of bubbles from 
fluid to ‘‘rigid’’ bodies as well as a criterion for this transition is most desirable, since such 
transition might be influenced by the magnitude of the pressure gradient. In the region of 
Reynolds numbers where the bubbles behave like rigid spheres, the pressure gradient probably 


has no effect on the drag coefficient. 


SUMMARY 


As the size of the bubbles was increased in the tests, a change in bubble shape from 
spherical to ellipsoidal to spherical cap shape was observed in all liquids. The volumes at 
which these transitions occur, however, varied with the properties of the liquid. For spheri- 
cal bubbles of given volume the results show that the viscosity of the liquid is the most im- 
portant property determining the rate of rise. For ellipsoidal bubbles, the surface tension 
assumes greater importance. Spherical cap bubbles rise independently of the properties of 
the liquid. 

The results show that the motion of air bubbles rising at their terminal velocity in a 
gravity field cannot be described completely by use of dimensionless parameters formed from 
the usual liquid properties (viscosity, surface tension, density), the equivalent radius of the 
bubble, the acceleration of gravity, and the terminal velocity. 

The drag coefficients of tiny spherical bubbles coincide with those of corresponding 
rigid spheres. With increase in bubble size, a decrease in the drag as compared to that of 
rigid spheres occurs in some liquids. This change in the drag is due to the development of 
circulation inside the bubble. The drag curves of the spherical bubbles rising in various 
liquids fall between two limiting curves, namely the drag curves of rigid and fluid spheres, 
respectively. It was not possible to determine a criterion for the transition region of the 
bubbles from ‘‘rigid’’ to fluid spheres. 


The region of ellipsoidal bubbles extends over different ranges of Reynolds numbers 
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for the various liquids. For liquids of low ‘‘@’’ number (say less than IO), a minimum in 
the drag curve is reached at Reynolds numbers of the order of 250. These minima occur near 
the transition from spherical to ellipsoidal shape. Such minima in the drag curve are not ob- 
tained for liquids of high ‘‘M’’ number. For the liquids used, transition to spherical caps is 
completed at a Weber number of about 20. 

The drag coefficients of spherical cap bubbles are independent of bubble size and 
have a constant value of 2.6. The rate of rise of these bubbles as a function of the equivalent 


radius is given by the experimentally determined relation: 


U=1.02Ygr, 


For bubbles (ranging in equivalent radius from 0.03 to 0.25 cm) rising in tap water, an 
_increased drag as compared to bubbles in clean (filtered or distilled) water was observed. 
The presence of certain surface-active substances in the water similarly increases the drag 
of bubbles (ranging in equivalent radius from 0.03 to 0.30 cm) as compared to bubbles in 
pure water. Beyond a certain critical concentration of these surface-active substances, an 
increase in concentration has relatively little influence on the drag of the bubbles. 

Tests to determine the effect of the container walls on the velocity of rise indicate 


the absence of such effect for the range of bubble volumes and container sizes tested. 
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APPENDIX 
RATE OF RISE OF GAS BUBBLES IN DISTILLED WATER 


It is of interest to compare the rate of rise of bubbles of different gases in distilled 
water. As pointed out under ‘‘Previous Experimental Work,’’ a number of tests have been 
conducted using oxygen bubbles. In Figure 30, the velocity of air and of oxygen bubbles in 
distilled water is given as a function of the equivalent radius of the bubble. The results 
indicate no significant difference in the rate of rise of the bubbles composed of oxygen and 


those of air. 
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Figure 30 - Terminal Velocity of Gas Bubbles in Distilled Water 
as a Function of Bubble Size 
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